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The relation between single particle and ensemble measurements is adressed for semiconductor 
CdSe nanocrystals. We record their fluorescence at the single-molecule level and analyse their emis- 
sion intermittency, which is governed by unusual random processes known as Levy statistics. We 
report the observation of statistical aging and ergodicity breaking, both related to the occurrence 
of Levy statistics. Our results show that the behaviour of ensemble quantities, such as the total 
fluorescence of an ensemble of nanocrystals, can differ from the time averaged individual quantities, 
and must be interpreted with care. 
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The relation between single particle and ensemble mea- 
surements is at the core of statistical physics, and be- 
comes crucial given that experiments are now able to re- 
solve individual nanometer-sized objects. This question 
is addressed here for semiconductor CdSe nanocrystals. 
The fluorescence properties of these colloidal quantum 
dots (QDs) have raised great attention due to their size- 
induced spectral tunability, high quantum yield and re- 
markable photostability at room temperature [Q, all of 
which make QDs a promising system for biological la- 
belling [g], single-photon sources [|| and nanolasers Q. 

When studied at the single molecule level, CdSe QDs 
share with a large variety of other fluorescent nanometer- 
sized systems j^, |], ffll the property of exhibiting fluores- 
cence intermittency ||^| . This means that the fluorescence 
intensity randomly switches from bright ('On') states to 
dark ('Off') states under continuous excitation. Although 
the very origin of the intermittency for CdSe QDs re- 
mains a matter of investigation, its statistical properties 
have been studied. For a given QD, the durations T on 
and r ff of the On and Off periods follow slowly decay- 
ing power-law distributions Pon(T on > r) = (ro/r) A ' on , 
P s(t b > t) = (70/7-/'°", where ^ on and are close 
to 0.5 Q [Uj 0. This behaviour extends over several 
orders of magnitude, from the detection integration time 
To up to hundreds of seconds, with very small dependence 
on temperature or excitation intensity. 

The crucial point for our analysis is that both /i on 
and //off are smaller than 1. In this case, the decay is 
so slow that the mean value of P on and P Q ff is formally 
infinite, and very long events tend to dominate the fluo- 
rescence signal, producing strong intermittency. The du- 
ration of the On and Off periods are thus governed by 
"Levy statistics", which have been encountered in vari- 
ous fields ||, [L§ P P, [l?], [ll], |(], such as laser 
cooling of atoms |14], dynamics of disordered |L5| and 



chaotic Jl7) systems, glassy dynamics or economics 
and finance |f20| . 

In this Letter, we show that single QD measurements 
can be used to explicitly compare ensemble- and time- 
averaged properties and explore some of the unusual phe- 
nomena induced by Levy statistics, such as statistical 
aging and ergodicity breaking. Using an epi-fluorescence 
microscopy set-up and a low-noise CCD camera, we si- 
multaneously recorded at room temperature the fluores- 
cence intensity of 215 individual QDs for duration of 10 
minutes with a time resolution of 100 ms (22). The blink- 
ing of the fluorescence intensity was observed for each 
QD detected in the field of the camera (Fig.|l|). Due to 
the binary behaviour of the blinking process, each in- 
tensity time trace was simply considered as a sequence 
of n On and Off times {t^ , t$ , t$ , t$ , . . . , ri™ ) , t$ } 
from which the distributions P on and P Q e were derived. 
In our measurements, the On and Off periods both fol- 
lowed power-law distributions After adjustment of 
the cumulative distributions of the On and Off periods 
for each of the 215 QDs, the exponents /i n and /K ff were 
estimated to be respectively 0.58 (0.17) and 0.48 (0.15), 
consistent with previous experiments pi] , [l2"| . For all pair 
of QDs, we also computed the Kolmogorov-Smirnov (KS) 
likelihood estimator |p4| to compare the On (resp. Off) 
distributions between each pair of QDs. For our set of 
data, the KS tests yield the same average value of 0.4 
(0.3) for both On and Off distributions, well above the 
value 0.05 usually considered as an inferior limit to as- 
sume that two datasets have identical distributions. In 
the following, the 215 QDs are therefore considered as 
statistically identical, with /i on =0.58 and /x o ff=0.48 |2^| . 

The first observation is that, for purely statistical rea- 
sons, the fluorescence of QDs is non-stationary, i.e. time 
translation invariance is broken in the intermittency pro- 
cess. This is best evidenced by studying the rate at which 
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Fig. 1: Fluorescence intermittency of a single CdSe nanocrys- 
tal measured over 10 minutes with a 100 ms time bins. Due 
to the broad distribution of the On and Off states, the signal 
is dominated by a few long events. 

the QDs jump back from the Off to the On state (a 
"switch on" event). For this purpose, we computed the 
ensemble average of the probability density s(9) to ob- 
serve a QD switching on between 9 and 8 + d9 after a 
time 9 spent in the Off state. For Off periods follow- 
ing a "narrow" distribution (with a finite mean value 
(Toff)), s(9) - also called the renewal density - would be 
independent of 9 and equal to l/(r ff). The situation is 
drastically changed for CdSe QDs : due to the fact that 
//off < 1, ( r off) no longer exists and s(9) is then ex- 
pected to decay as 6>~( 1_AIoff ) |Q. This means that as 
time grows, the switch on events occur less and less fre- 
quently. Fig.||a shows that our data match these theoreti- 
cal predictions : the measured value of s(9) decreases like 
9~ a , with a = 0.5 in agreement with the value (1 — ^ ff) 
expected from our measurement of /^ ff - 

This non-stationary behaviour can be understood by 
considering, for each QD, the quantity : 

N 
i=l 

i.e. the total time spent in the Off state during the N 
first Off periods (FigJ|b) . Assuming that the are in- 
dependent and T g having no mean value, the sum of N 
such independent random variables must be evaluated by 
means of the Generalized Central Limit Theorem (see e.g. 
|0). This theorem states that 9(N), instead of scaling 
as N, grows more rapidly, as N 1 ^ " . As shown in Fig.^D, 
the sum 9(N) is dominated by few events. This central 
property, distinctive of Levy statistics, means that, as 
time grows, one observes long events that are of the or- 
der of 9(N) itself HJ. Hence, the probability for a QD 
to switch on decreases with time : the system ages 
and the signal is non-stationary. 

To test the assumption that the Off events are inde- 
pendent and to gain further insight into this aging effect, 
we computed the persistence probability Hq(9,0 + 9'), 
defined as the probability that no switch on event oc- 
curs between 9 and 9 + & . In the case of independent 
Off periods with an exponential distribution (with mean 




Fig. 2: Statistical aging measured from a sample of 215 QDs. 
(a) Logarithmic plot of the probability density s(9) for a QD 
to jump in the On state after having spent a total time 9 
in the Off state. The solid line is a power-law adjustment 
9~ a with a — 0.5. (b) Evolution of the total time spent in 
the Off state 9(N) with the number N of Off periods for a 
given QD ; the sum 9(N) is dominated by a few events of the 
order of 9(N) itself, (c) Persistence probability U o (0,0 + 9') 
measured from the set of 215 QDs for 6 = 0.1 s (A), 6 = 0.5 
s (o), 9 = 1 s (x), 9 = 2 s (□), 9 = 4 s (*), (9 = 8 s (a). 
Uo(0, 9+9') depends on 9, indicating that the process is aging, 
(d) Persistence probability Ho(9, 9 + 9') measured for 9 and 9' 
ranging between and 10 s with 0.1 s time bins and expressed 
as a function of 9/(9 + 9') (+). Each point corresponds to the 
average over 200 adjacent data points. The solid lines are the 
theoretical predictions for exponents /i ff = 0.5, /u g = 0.55 
and jUog = 0.6. 

value (T ff)), Ho(9, 9 + 9') is independent of 9, and given 
by e~ e '( Tolf ', illustrating that the switching process is 
stationary. The computation of Ho from our dataset re- 
veals a completely different pattern : the probability that 
no switch on event occurs within a given duration 9' de- 
creases with 9 (Fig.||c), consistent with the behaviour 
of s(9). Furthermore, Hq(9,9 + 9') is found to depend 
only on the reduced variable 9/ (9 + 9'), and to vanish for 
9/(9 + 9') close to (Fig^d). This result proves that one 
has to wait a time 9' of the order of 9 to have a chance 
to observe a switch on event, in qualitative agreement 
with the fact that the largest term of the sum 9(N) is of 
the order of 9(N) itself. Quantitatively, for independent 
Off events distributed according to a Levy distribution 
Poff with exponent /j, s , the persistence probability is 
expected to read : 

,0/(0+0') 

U (9,9 + 9')= / l3^ ttA . M (u)du 
Jo 

where (3 is the beta distribution on [0,1] jl9[ Our 
data follow this prediction with fi g = 0.55, in agreement 
both with our previous estimations of /z ff (Fig.||d) and 
with the assumption that the Off events are independent. 
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These results show that the aging effect has a pure sta- 
tistical origin and is not related to an irreversible process 
(such as photo-destruction). Due to the statistical prop- 
erties of Levy distributions, non-stationarity emerges de- 
spite the time-independence of the laws governing the 
microscopic fluorescence process. 

From a more general standpoint, this non-stationary 
behaviour has also profound consequences on basic data 
interpretation, such as the ensemble-averaged total fluo- 
rescence emitted by a population of QDs. We illustrated 
this by studying & on (t), the fraction of QDs in the On 
state at a given time t (Fig.||a). In the context of Levy 
statistics, the time evolution of $ on (i) is intimately linked 
to the relative amount of time spent in the On and Off 
states for each QD. Qualitatively, the Off events tend to 
be dominant whenever /i, g < fi g since 9(N) = T oS 
grows faster than its counterpart 6(N) = J2iLi T on • 
When analysed in a more quantitative way, the frac- 
tion $ on (t) can be shown to decrease asymptotically as 
f**o«-j*on [Q. Experimental results confirm this analysis : 
3\>n(i) decays as t~@, with an exponent [3 = 0.13 indeed 
consistent with the previous determination of /j, on and 
pL g (Fig.||a). We also observed that the average signal 
over the whole CCD detector - i.e. the sum of the fluo- 
rescence of all the QDs - decays as i~ 018 , in agreement 
(within experimental uncertainty) with the fact that time 
increasing, less and less QDs are in the On state, caus- 
ing the total fluorescence to decrease like $ O n(0 (Fig.||b). 
Importantly, we also observed that this fluorescence de- 
cay is laser-induced and reversible : after a continuous 
laser illumination of 10 minutes, leaving the sample in 
the dark for about 10-15 minutes systematically lead to 
a complete recovery of its initial fluorescence. This con- 
firms that this decay is again purely statistical, and not 
related to an irreversible bleaching of the QDs. 

Our final observation focuses on non-ergodic aspects of 
random processes driven by Levy statistics. Single par- 
ticle measurements allow one to compare directly < & n(i) 
and the fraction of time $on(0 i— ► t) spent in the On 
state between and t for the zth QD. This provides a di- 
rect test of the ergodicity of the QD fluorescence. While 
the ensemble average <& n(i) decays deterministically as 
t -013 (Fig.pk), each time average widely fluctuates over 

(i) 

time and for a given t, the values of <&on are broadly 
distributed between and 1, even after a long time of 
integration (Fig.||c) . To study the behaviour of time av- 
erages, we calculated the relative dispersion oy(t) of the 
time averages at time t, where ay(t) corresponds to the 
standard deviation of the distribution of $on (0 i— » t) over 
the set of QDs, divided by its mean value. Fig.||d shows 
that a r (t) does not decay to zero, and is still of order 1 
on the experimental time scale. Therefore, even for long 
acquisition times, the fluctuations of the time averages 
from QD to QD remain of the order of the time averages 
themselves and do not vanish as expected for ergodic sys- 
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Fig. 3: Non-stationarity and non-ergodicity in a sample of 
QDs. (a) Time evolution of the fraction $ on (i) of QDs in the 
On state at time t (V). $ on (t) decays as t~° 13 = 
(solid line), (b) Time evolution of the total fluorescence sig- 
nal emitted by the sample (□) : the darkening effect follows a 
t -0 ' 18 power-law decay (solid line), (c) Typical time evolution 
of $on(0 i— ► t) - the fraction of time spent in the On state 
between and t - for 7 QDs. The time averages are widely 
fluctuating, even in the long integration time limit, (d) Evo- 
lution of the relative dispersion o r {t) of $in (0 i— ► t) at time t 
over the ensemble of QDs (o). As time grows, a r (t) tends to a 
constant value, illustrating that the time averages trajectories 
do not converge to any asymptotic value. 

terns. These data indicate ergodicity breaking : due to 
rare events with a duration comparable to the total ac- 
quisition time, there is no characteristic time-scale over 
which physical observables can be time-averaged. Even 
for long acquisition time, $on(0 i— > t) does not converge 
and no information on the ensemble value $ n can be 
obtained by time averaging an individual trajectory. 

While we found that accurate estimates of fi on and ^t ff 
are essential to analyze and predict the statistical prop- 
erties of the fluorescence, the microscopic origin of these 
broad distributions is not yet established. Distributions 
of Off times are sometimes attributed to distributions of 
static traps from which the charge of an ionised QD es- 
capes by tunnelling effect Jll], In these models, the 
value of ^off strongly depends on microscopic characteris- 
tics of the QDs, and it is not clear how this is compatible 
with the statistical homogeneity of the different QDs sug- 
gested by the KS test. The dynamic changes of the parti- 
cle environment are also often invoked to account for the 
fluctuating emission of the QD [ |i~2] , [2S| |. Some authors 
have thus suggested models in which the trap for the 
charge of the ionised QD follows a random walk in a ID 
parameter space, yielding a universal value 1/2 for jj, s 
Jj~2[ . However both of these models (static and dynamic) 
have yet to be more thoroughly tested. Since intermit- 
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tency is an ubiquitous process at the nanometer scale, 
some of the arguments discussed here for QDs might also 
apply to other systems. In particular, our analysis shows 
that non-stationary behaviour of the fluorescence - some- 
times attributed to photochemical processes - can also 
have purely statistical origins (such as statistical aging). 
Recent evidence have shown that this may be the case in 
a system as microscopically different from QDs as green 
fluorescent proteins |3(]]. In this respect, aging and non- 
ergodicity might be an important pattern when studying 
single nanometer-sized objects in complex environments. 

In conclusion, our experimental results show that en- 
semble averaged fluorescence properties of individual 
CdSe QDs are deeply affected by the non-standard sta- 
tistical properties of the Levy statistics governing the 
blinking process. We found that a population of QDs ex- 
hibit statistical aging. Hence, despite the blinking statis- 
tics are time- independent, the fluorescence emitted by an 
ensemble of QDs under continuous laser excitation is non- 
stationary. Our data also evidence that due to the scaling 
properties of Levy statistics, CdSe QDs are non-ergodic 
systems : time- and ensemble-averaged properties do not 
coincide anymore, in full contrast with usual assumptions 
when studying nanoscale emitters. 
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